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Abstract. Asteroseismology has been recognized for a long time as a very powerful mean to 
probe stellar interiors. The oscillations frequencies are closely related to stellar internal structure 
properties via the density and the sound speed profiles. Since these properties are in turn tightly 
linked with the mass and evolutionary state, we can expect to determine the age and mass of a 
star from the comparison of its oscillation spectrum with the predictions of stellar models. Such 
a comparison will of course suffer both from the problems we face when modeling a particular 
star (for instance the uncertainties on its global parameters and chemical composition) and from 
our general misunderstanding of the physical processes at work in stellar interiors (for instance 
the various transport processes that may lead to core mixing and affect the ages predicted 
by models). However for stars where observations have provided very precise and numerous 
oscillation frequencies together with accurate global parameters and additional information (as 
the radius or the mass of the star if it is member of a binary system, the radius if it observable 
in interferometry or the mean density if the star is an exoplanet host), we can also expect to 
better constrain the physical description of the stellar structure and transport processes and to 
finally get a more reliable age estimation. 

After a brief survey of stellar pulsations, we present some general seismic diagnostics that can 
be used to infer the age of a pulsating star as well as their limitations. We then illustrate the abil- 
ity of asteroseismology to scrutinize stellar interiors on the basis of a few exemples. In the years 
to come, extended very precise asteroseismic observations are expected, either in photometry or 
in spectroscopy, from present and future ground-based (HARPS, CORALIE, ELODIE, UVES, 
UCLES, SIAMOIS, SONG) or spatial devices (MOST, CoRoT, WIRE, Kepler, PLATO). This 
will considerably enlarge the sample of stars eligible to asteroseismic age determination and 
should allow to estimate the age of individual stars with a 10-20% accuracy. 
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1. Introduction 

As can be seen in Fig.[TJ pulsating stars are presently observed in nearly each re- 
gion of the HR diagram where stars are observed. Stellar pulsations may be excited by 
different mechanisms, in a large range of amplitudes. Self-excitation results from the 
K-mechanism which drives pulsations either in the Hell ionisation zone (oscillations in 
S Scuti, Cepheids, RR Lyrae or DB white dwarfs), in the HI and Hel ionisation zones 
(oscillations in roAp, Miras and irregular variables) or in the metal opacity bump (see 
e.g, Dziembowski et al. 1993 ; Pamiatnykh 1999 for oscillations in (3 Cephei and slowly 
pulsating -SPB- stars). In 7 Doradus stars, excitations are interpreted as due to the 
convective blocking of the radiative flux (see e.g, Guzik et al. 2000 ; Warner et al. 2003 ; 
Dupret et al. 2004). Finally, in stars with significant convective envelopes, the so-called 
solar-like oscillations are interpreted as resulting from stochastic excitation by turbulent 
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Figure 1. Pulsators in the HR diagram (courtesy J. Christensen-Dalsgaard). 

convective motions. This can occur on the main sequence (Christensen-Dalsgaard 1982) 
and in subgiants and red giants (Dzicmbowski et al. 2001). 

Oscillations are observed either by spectroscopic measurements (velocity variations) or 
by photometric ones (intensity variations). While only one or two modes but with large 
amplitudes are observed in Cepheids, W Virginis or RR Lyrae stars, a large number 
of modes with small amplitudes are observable in solar-like pulsators, S Scuti, roAp, [5 
Ccphci, SPB, 7 Doradus and white dwarfs pulsators. In solar- like pulsators the variations 
in velocity amplitudes range from less than 1.0 to about 50 m.s -1 (Sun: 0.2 m.s -1 ) while 
relative intensity changes vary from 10~ 6 to 10~ 3 (Sun: 4. 10~ 6 ) and the amplitude 
scales as (^) where a is in the range 0.7 — 1.0 (Kjeldsen and Bedding 1995 ; Samadi 
et al. 2005). Amplitudes in 8 Scuti, roAp, (3 Cephei, SPB, 7 Doradus and white dwarfs 
pulsators range from 10 -3 to a few 10 _1 mag. 

2. Characteristics of stellar oscillations 

Each oscillation cigenmode of frequency v n ,£,m consists in a radial part characterized 
by the radial order n (n is the number of nodes along the stellar radius) and in a surface 
pattern characterized by a spherical harmonic Ye^ m where I is the mode angular degree 
(£ is a measure of the number of wavelengths along the stellar circumference) and m is 
the azimutal order (number of nodes along the equator). 

In the resolution of the equations of stellar oscillations two characteristic frequencies 
appear, in addition to the cut-off frequency. First, the acoustic (Lamb) frequency Si — 
[£(£ + l)]2(c/r) which is a measure of the compressibility of the medium and depends on 
the adiabatic sound speed c = {T\P/p)2 whereTi = (91nP/91np) a( j is the first adiabatic 
coefficient. Second, the buoyancy or Brunt- Vaisala, frequency Nbv, such that iV| v = 

9 ( fT dl d r P ~ ) • ^bv i s positive in a radiative region where it corresponds to the 




Figure 2. Propagation diagrams for a solar model (left), a 1.4 Mq model on the subgiant branch 
(centre) and a 2 M© model (5 Scuti) at the end of the MS (right). The frequency domains of 
solar-like oscillations (high degree p-modes) is located in-between the dotted lines. 



oscillation frequency of a perturbed fluid element while it is negative in convective regions. 
It is worth noticing that for ideal gases, c oc (T/ p) 2 and 7V| V ~ (g 2 p / 'P)(V ad — V + V M ) 
where p is the mean molecular weight and V M is its gradient. 

It can be shown that different waves can propagate inside a star: modes with frequencies 
higher than both Si and iVav correspond to standing sound waves (acoustic pressure 
modes or p-modes) while modes with frequencies lower than both Sg and -/Vbv correspond 
to standing gravity waves (g- modes). Otherwise modes are evanescent. In absence of 
rotation, the frequency of an oscillation mode only depends on two parameters (density 
p and Ti or equivalently p and sound speed). Therefore frequencies change with mass, 
evolution (age) and chemical composition. 

In Fig. [2 we show propagation diagrams. The left figure shows that in the present 
solar interior the domains of propagation of p and g-modes are well separated. The 
central figure shows a model of a star of 1.4 Mq on the subgiant branch. In this star a 
/i-gradient has been built in the radiative core during the main sequence (MS) evolution 
and the core has become more and more dense on the subgiant branch. As a result, the 
Brunt- Vaisala frequency shows a central peak; this allows g-modes to propagate in a 
domain of frequency corresponding to the p-modes domain. Due to the closeness of the 
propagation regions (see the £ = 1 mode on the figure), a g-mode can interact with a 
p-mode of the same frequency through the so-called avoided crossing and we therefore 
expect to observe p-g mixed modes in this kind of stars (Aizenman et al. 1977). This 
situation is also expected to happen in S Scuti stars (right figure), on the MS, because of 
the building of a /x-gradient in the radiative regions just above the receding convective 
core. More generally, depending on the mass and evolutionary state, different kinds of 
modes are predicted to be observable: for instance in (3 Cephei and 5 Scuti stars low 
order p and g modes are expected with periods in the range 2-8 hours for the former and 
30 min-6 hours for the latter, in SPB stars and 7 Doradus stars high order g-modes are 
expected with periods in the range 15 hours-5 days for the former and 8 hours-3 days 
for the latter while in solar-like pulsators high order p-modes are expected with periods 
in the range a few minutes (MS stars) to a few hours (red giants). 



3. Age diagnostics through p-modes 

In the case of the Sun, more than 10 5 modes can be observed, of all degrees, with an 
excellent precision. It is therefore possible to derive the density, the sound speed or the 
internal rotation rate by inversion of the frequencies which in turn gives access to very 
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Figure 3. Detail of the solar oscillations amplitude spectrum, as observed by VIRGO on SOHO 
(Figure from Bedding & Kjeldsen 2003). Each frequency peak is labelled with the value of the 
mode order and degree and the large and small separations are indicated. 

precise information on the solar properties such as the age, helium content and depth 
of the outer convective envelope (see e.g. Christensen-Dalsgaard, these proceedings). On 
the other hand, in the case of solar-like pulsators, only a few modes of low degree (mainly 
I = 0, 1, 2, 3) are accessible to observations. However, as discussed below, this can provide 
valuable information on the stellar interior. 

3.1. Asymptotic theory: frequency combinations and asteroseismic diagrams 

Figure [3] shows a detail of the solar power spectrum where the modes have been identified 
by their order and degree. The regular pattern seen is indeed predicted by the asymptotic 
theory of stellar oscillations (Vandakurov 1967, Tassoul 1980) from which the frequency 
of a mode of high order n and of degree I <C n is given by v n _g = Av (n + \l + e n A — 
£(£ + 1)Dq where e rit i is a function sensitive to surface physics, but only weakly sensitive 
to the order and degree of the mode, and Av and Dq will be explicited below. As a 
consequence the difference in frequency between two modes of consecutive orders and 
same degree is approximately constant and given by Ai> ~ v n +\,i — v n j = Avi while the 
difference in frequency between two modes of consecutive orders and degrees differing by 
two units is Svi = v ai £ — v n -i.e+2 = &(£ + 6)Dq. The differences Av and Sv are called the 
large and small frequency separations, respectively. 

It can be shown that the large frequency separation can be approximated by Av ~ 
(2 §q ^r) _1 - Therefore it is a measure of the inverse of the time it takes a sound wave 
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Figure 4. Left: the change of the mean density with evolution on the MS for 3 values of the 
stellar mass. Centre and right: the sound speed profile (centre) and the dc/dr profile (right) in 
1 M Q models, on the ZAMS, in the middle of the MS and at the end of the MS (TAMS) 
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to travel across the star and it scales as (M/i? 3 ) 1 / 2 , that is as the mean density. As 
shown in Fig. [4] left , for stars on the MS, the mean density increases as mass decreases 
and at given mass it decreases as evolution proceeds; we therefore expect the large 
separation to be quite sensitive to the mass and slightly to the evolutionary state on 
the MS. On the other hand, the small frequency separation depends on the value of 
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and is therefore sensitive to the integral over the star of 
the sound speed gradient weighted by 1/r. During MS evolution, both the temperature 
T and the mean molecular weight [i increase in stellar cores but the relative increase in 
\i exceeds that in T so that the net effect is a decrease of the sound speed towards the 
center (see Figs. |4] centre and right). This in turns leads to a decrease in Do and therefore 
of the small separation as evolution proceeds making Sv a valuable age estimator. 

Fig. [5] left recalls how the age of a star is poorly constrained in those regions of the HR 
diagram where the isochrones are degenerate (for instance in the regions close to the zero 
age main sequence of low mass stars and in the turnoff regions) . On the contrary if we plot 
in a diagram the large separation versus the small one (Christensen-Dalsgaard 1988), or 
rather the ratio — 5vi / Avi of the small to the large separation which is less sensitive 

to surface effects (Roxburgh & Vorontsov 2003), we find that the degeneracy is removed 
(see Fig. [5] right). Therefore a precise measure of p-modes frequencies should provide a 
valuable diagnostic of the age (and mass) of solar-like oscillators. Indeed Kjedsen et al. 




Figure 5. Left: the inversion of isochrones in the HR diagram doesn't allow to estimate the age 
of non evolved stars (the age of a G-star similar to the Sun is found to be in the range 0.5-10 
Gyr!) while it gives poor accuracy on the age of stars close to the turnoff (the age of a typical 
A-star lying at the turnoff is found to be in the range 0.6-1 Gyr). Large ellipses correspond 
to the present uncertainty on the HR diagram position while small ones correspond to the 
accuracy achievable at GAIA time (see e.g. Lebreton 2005). Isochrones have been obtained 
with the CESAM code (Morel & Lebreton 2008). Right: Asteroseismic diagram in which the 
abscissa is the mean large frequency separation Av of £ — modes of different orders and the 
ordinate is the mean ratio of the small separation (between 1 = 1 and I = modes) to the large 
separation. The frequencies have been calculated for standard models of different masses and 
evolution stages on the MS. Models have been obtained with either the CESAM or the CLES 
code (Morel & Lebreton 2008, Scuflaire et al. 2008a) and frequencies with the LOSC oscillation 
code (Scuflaire et al. 2008b). The frequency means have been estimated following Mazumdar 
(2005). The error bar located at the solar position corresponds to the one expected from CoRoT 
or Kepler missions measurements, i.e. to a relative accuracy on the individual frequency of 10 -4 . 
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(2008) have estimated that with an accuracy on observed frequencies of a few 0.1 /jRz 
(typical accuracy of the CoRoT and future Kepler missions), the age of a solar type 
oscillator can be determined to better than « 10% of its MS lifetime. 

In the age estimation it is important to evaluate the uncertainties that result either 
from the uncertainty on the determination of the stellar fundamental parameters (as the 
chemical composition) or from the weak knowledge of several aspects of the physics enter- 
ing stellar models calculation. For instance the description of convection, the occurence 
of either microscopic or turbulent diffusion of chemical elements or the estimate of the 
size of the mixed core -which depends on processes like overshooting or rotational mixing- 
will all affect the calculated frequencies. We have estimated how the age estimated from 
an asteroseismic diagram ((7-0,2), (A^o)) changes with stellar model inputs. We find that 
changing the mixing-length parameter from the value ojmlt = 1-65 to 1.80 changes the 
age at a given evolutionary stage on the MS (same central hydrogen content X c ) by less 
than 1.5% while introducing microscopic diffusion in the model calculation changes the 
age by less than ss 5%. Changing the initial helium abundance in mass fraction from 
Y = 0.27 to 0.26 modifies the age at constant X c by less than w 10%. On the other 
hand, as illustrated in Fig. [5] (left) changing the initial metallicity in mass fraction from 
the value Z = 0.02 to 0.01 changes the age by 15 to 30% while changing the overshooting 
parameter from a ov — 0.0 to 0.2 has huge effects on the age estimate for stars possessing 
convective cores in the late part of their MS evolution (see Fig. [6] right). 

To improve the determination of the age based on asteroseismic diagrams, it is therefore 
crucial to improve the observational determination of their metal content (this requires 
high resolution spectroscopy and further improvements in model atmospheres). Also it 
is very interesting to get additional (preferably independent) information on their fun- 
damental parameters. This is achievable for stars observable in interferometry (direct 
access to the radius), for members of binary systems (access to mass and/or radius), or 
for hosts of exoplanets (access to the mean density from the planetary transit). Finally, 
as discussed in the following, the full interpretation of seismic data may allow to get 
further information, such as the surface helium abundance or the size of the mixed core. 

3.2. Deviations from the asymptotic theory: signature of sharp features in the interior 

Rapid variations of physical quantities in the interior of a star can be detected by seismic 
analysis. Indeed, any sharp feature inside the star will produce an oscillatory signal in 
the frequencies of the form bv ~ v4cos(2r m ix) + (j)) where r m = is the acoustic 

depth of the abrupt feature (see e.g. Gough 1990, Monteiro et al. 2000). 

For instance, the second He ionisation in the outer convection zone of solar-like stars 
produces a depression in the profile of the adiabatic exponent Ti. The higher the helium 
abundance in the convection zone, the deeper the Ti depression. The sound speed profile 
locally carries the signature of the Ti depression which results in an oscillatory signal 
in the frequencies as a function of radial order. Basu et al. (2004) have shown that this 
signal, that can be extracted from low degree p modes frequencies, can be used to infer the 
helium abundance in the envelope of solar-like pulsators. They find that, with a relative 
accuracy of 10~ 4 on the frequencies (typical of C0R0T or Kepler missions), it should be 
possible to infer the envelope helium abundance (Y in mass fraction) of pulsators of mass 
in the range 0.8 — 1.4 M@ with an accuracy of AY « 0.02 provided their radius or mass 
can be estimated independently. Also, it should be possible to infer the position of the 
basis of the convective envelope (see e.g. Ballot et al. 2004) or the boundary of the mixed 
core (see e.g. Mazumdar et al. 2006) in a similar way. The latter depends on the extra 
mixing processes that can take place at the border of the convective core (overshooting, 
rotationally induced mixing) and it is a crucial data for age estimation. 
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Figure 6. Left: effects on the asteroseismic diagram of a change of chemical composition. In 
black models, with Y = 0.27, Z = 0.02; in dark grey, models with Y = 0.26, Z = 0.02 and 
in light grey, models with Y = 0.27, Z — 0.01. Right: Effect of overshooting with in black the 
models without overshooting and, in grey, models calculated with a ov = 0.2. 
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Figure 7. Left: Hydrogen abundance in the core of GMq models on the main sequence (upper 
panel) at X c ~ 0.5 (dotted line) and at X c ~ 0.3 (dashed line). The convective core recedes 
during the evolution leaving behind a chemical composition gradient and, a sharp feature in the 
Brunt- Vaisala frequency (lower panel). Right: Period spacing for the same models as a function 
of radial order k. Horizontal dotted lines represent constant period spacing as predicted by the 
asymptotic approximation, i.e. AP ki < oc n ( ) (Miglio et al. 2008). 

4. Age diagnostics through g-modes and p-g mixed modes 

The properties of g-modes are determined by the behavior of the Brunt- Vaisala fre- 
quency in the stellar interior. These modes are sensitive probes of the chemical gradients 
resulting from the combined effect of nuclear burning and convective mixing, and hence 
they are potential diagnostics of the stellar evolutionary stage (and therefore of age) . 

In the case of high order low degree g-modes, such as those expected in SPB (4 — 
7 Mq) and 7 Doradus stars (1.4 — 1.8 M Q ), the first order asymptotic theory predicts 
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Figure 8. Diagram in which the abscissa 
is the main period spacing for I = 1, and 
the ordinate is the period of the oscillatory 
behavior of the period spacing (see Fig. [Jj 
right) in terms of the difference in radial or- 
der between two consecutive trapped modes. 
Dotted lines represent the variation of these 
quantities along the MS evolution for mod- 
els with masses between 1.4 and 8 Mq, initial 
chemical composition Xo — 0.70, Zo = 0.02, 
and overshooting parameter «ov = 0.0. 
Solid lines (grey) connect models at the same 
evolutionary state, i.e. with the same value 
of the central hydrogen mass fraction. 
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that the periods of the modes should be equally spaced (Tassoul 1980) with a period 
spacing for two modes of consecutive order fc and k + 1 and same degree £ given by 

APk,^ = Pk+i,i ~~ Pk,i — t - Hoi with the buoyancy radius Hq = (^J^ ^ N f v ^ drj . This 
approximation is, however, no longer valid when variations of Nbv on a length scale 
smaller than the oscillation wavelength are present. That may occur in MS stars due 
to the building of a /Lt-gradient at the outer border of the receding convective core. The 
corresponding sharp feature in TVbv can lead to a resonant condition such that modes of 
different order k are periodically confined in the region. These modes have a different 
period and as a consequence the period spacing presents an oscillatory behavior rather 
than the constant value that would be expected in a model without sharp variations 
in iVev (Fig- E]). By using a second order approximation, it is possible to derive an 
analytical expression relating the amplitude and periodicity of the oscillatory component 
to the location and sharpness of the feature in A^bv (Miglio et al. 2008 and references 
therein). This periodicity expressed in terms of radial order is given by the ratio between 
the buoyancy radius of the star, and the buoyancy radius of the V M region: Afc = H^/TIq, 
and corresponds to the difference between radial orders of two consecutive trapped modes 
in the V M region. As shown in Figs. [7J and [5J the value of Afc is very sensitive to the 
evolutionary state of the star: for non evolved models (X c ~ Xq) Afc — » oo and the 
period spacing is almost constant as predicted by the asymptotic theory but when the 
stars evolve the periodicity of the oscillatory component in AP decreases. 

(3 Cephei and S Scuti show p and g-modes of low radial order and, for stars evolved 
enough, also p-g mixed modes (see. Sect. HJ. In this case, the quantitative predictions 
of the second order asymptotic approximation cannot be used, but they are still able to 
qualitatively describe the properties of g-modes. The periods of g-modes depend on the 
location and shape of the chemical composition gradient and therefore so do the frequen- 
cies of mixed modes. Since the frequency separation between consecutive modes varies 
rapidly during the avoided crossing phenomenon, the detection of mixed modes together 
with that of pure p-modes provides an important constraint to the stellar evolutionary 
state (see for instance Pamyatnykh et al 2004). 

Finally for low mass stars that leave the MS and are evolving as sub-giants, the in- 
creasing central condensation together with the chemical composition gradient lead to a 
large increase of A^bv- As a consequence, mixed modes may appear in the frequency do- 
main of solar-like oscillations. Again, the frequencies of the modes undergoing an avoided 
crossing are direct probes of the evolutionary status: models on the MS do not present 
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mixed modes and the large frequency separation is almost constant (as predicted by the 
asymptotic theory). On the contrary, as more evolved models are considered, the avoided 
crossing effects become more important and the large frequency spacing for non-radial 
modes becomes more irregular (see e.g. Di Mauro et al. 2004, Miglio et al. 2007). 

In all the cases considered above, the frequencies of high-order g-modes and g-p mixed 
modes are sensitive to the location and shape of the chemical composition gradient. 
Therefore, these modes are sensitive probes of the evolutionary state and also of the 
properties of the transport processes acting in stellar cores, for instance, convective over- 
shooting or diffusive mixing. Nevertheless, to clearly discriminate between different sce- 
narios, high accuracy of the classical observables such as effective temperature, luminosity 
and chemical composition (and masses for binary stars) would also be required. 

Age estimates of stars in the galactic disc are crucial for galactic evolution studies. 
Considering the progress expected on the determination of both the seismic and classical 
observables, Lebreton et al. (1995) (see also Lebreton 2005) find that we can reasonably 
expect to reduce the uncertainty on the age determination of disc A and F stars which is 
presently of about 35 — 40% (of which 13 — 24% result from the uncertainty on the size 
of the mixed cores) to a level of ~ 10 — 15%. 

5. Prospects and Conclusions 

We are presently in a very fruitful period where asteroseismic data are being obtained 
both from the ground and from space to an increasing level of precision for growing 
samples of stars in wide ranges of masses and evolutionary stages. The CoRoT satellite 
has observed for more than one year now. It will observe a total of more than 100 
targets (solar-like stars, 7 Doradus, <5 Scuti, j3 Cephei, giants) in the so-called Seismo 
field and for stars observed during Long Runs of 150 days, a frequency accuracy as 
good as a few 10~ 7 Hz will be reached (see Michel et al. 2006). First C0R0T observations 
are currently under analysis and interpretation. Up to now, three solar-like pulsators 
have been observed (see e.g. Michel et al. 2008a) and for the first one HD49933, observed 
during the first 60 days Intermediate Run and then during a 150 days Long Run, the large 
frequency separation has been measured and a first frequency identification for modes of 
degrees £ = 0, 1, 2 has been proposed (Appourchaux et al. 2008). Further very interesting 
results are being presently obtained related for instance to solar-like oscillations in giant 
F and G stars and to the very rich oscillation spectra of S Scuti stars just revealed by 
C0R0T (see e.g. Michel et al. 2008b). 

The Kepler photometry space mission to be launched in March 2009 will provide 
seismic data, at the same level of precision than C0R0T, for a large number of stars, hosts 
of exoplanets. The information on the stellar mean density, coming from the observation 
of the exoplanet transit, will add constraints for the age determination. Asteroseismic 
observations are presently obtained in spectroscopy with HARPS, CORALIE, ELODIE, 
UVES, UCLES and in photometry with MOST, C0R0T, WIRE while other projects are 
about to begin or are currently under study either in spectroscopy (SONG, SIAMOIS) or 
in photometry (Kepler, PLATO). 

In parallel, much progress will be done in the determination of the stellar global pa- 
rameters. For instance GAIA (ESA 2000; Perryman et al. 2001), to be launched in 2011, 
will make astrometric measurements, at the micro-arc second level together with photo- 
metric and spectroscopic observations of a huge number of stars covering the whole range 
of stellar masses, compositions and evolution stages. Important progress is also expected 
from interferometric and high resolution spectroscopy measurements (for instance with 
VLT-VLTI, CHARA, KECK, JWST etc.). 
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Asteroseismology provides diagnostics of stellar evolution all across the HR diagram. If 
we consider both the present and expected future accuracies on the observed frequencies 
and the future improvements expected in the determination of the stellar global param- 
eters, it seems reasonable to forecast that asteroseismology will allow to determine the 
ages of individual stars with a 10 — 20% accuracy. This will certainly lead to significant 
progress in the understanding of the history and evolution of galaxies. 
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